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We study the influence of entanglement and correlated noise using cor- 
^ . related amplitude damping, depolarizing and phase damping channels on 

' the quantum Stackelberg duopoly. Our investigations show that under 

the action of amplitude damping channel a critical point exists for unen- 
tangled initial state as well, at which firms get equal payoffs. The game 
becomes a follower advantage game when the channel is highly decohered. 
Two critical points corresponding to two values of the entanglement angle 
are found in the presence of correlated noise. Within the range of these 
' limits of entanglement angle, the game is follower advantage game. In case 

, of depolarizing channel, the payoffs of the two firms are strongly infiuenced 

' by the memory parameter. The presence of quantum memory ensures the 

existence of Nash equilibrium for the entire range of decoherence and en- 
^— «j . tanglement parameters for both the channels. A local maximum in the 

' payoffs is observed which vanishes as the channel correlation increases. 

, Moreover, under the influence of depolarizing channel, the game is always 

a leader advantage game. Furthermore, it is seen that phase damping 
channel does not effect the outcome of the game. 
^ ; PACS: 02.50.Le, 03.67.-a, 05.30.-d. 

\^ • Keywords: Quantum channels; correlated noise; Stackelberg duopoly. 
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1 Introduction 

Game theory is the mathematical study of interaction among independent, self 
interested agents. It emerged from the work of Von Neumann [1], and is now 
used in various disciplines like economics, biology, medical sciences, social sci- 
ences and physics |21 13] . Due to dramatic development in quantum information 
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theory [?], the game theorists [5-12] have made strenuous efforts to extend the 
classical game theory into the quantum domain. The first attempt in this direc- 
tion was made by Meyer |13| by quantizing a simple coin tossing game. Applica- 
tions of quantum games are reviewed by several authors [14-18]. A formulation 
of quantum game theory based on the Schmidt decomposition is presented by 
Ichikawa et al. [T^. Recently, Xia et al. (301 [21] have investigated the quantum 
Stackelberg duopoly game under the influence of decoherence and have found 
a critical point for the maximally entangled initial state against the damping 
parameter for the amplitude damping environment under certain conditions. 

In practice no system can be fully isolated from its environment. The in- 
teraction between system and environment leads to the destruction of quantum 
coherence of the system. It produces an inevitable noise and results in the loss 
of information encoded in the system |22| . This gives rise to the phenomenon 
of decoherence. Quantum information are encoded in qubits during its trans- 
mission from one party to another and require communication channels. In a 
realistic situation, the qubits have a nontrivial dynamics during transmission 
because of their interaction with the environment. Therefore, a party may re- 
ceive a set of distorted qubits because of the disturbing action of the channel. 
Studies on quantum channels have attracted a lot of attention in recent past 
P51 [21] • Early work in this direction was devoted mainly to memory less channels 
for which consecutive signal transmissions through the channel were not corre- 
lated. In the correlated channels (i.e. the channels with memory), the noise 
acts on consecutive uses of the channel. The effect of decoherence and corre- 
lated noise in quantum games have produced interesting results and is studied 
by number of authors [23-28] . 

In this paper, we study the effect of correlated noise introduced through 
amplitude damping, phase damping and depolarizing channels parameterized 
by the decoherence parameters pi and p2 and the memory parameters /x^ and 
/L<2j on the quantum Stackelberg duopoly game. The decoherence parameters pi 
and the memory parameters /Uj range from to 1. The lower and upper limits 
of decoherence parameter pi correspond to the undecohered and fully decohered 
cases, respectively. Whereas the lower and upper limits of memory parameter 
/Zj correspond to the uncorrelated and fully correlated cases, respectively. It is 
seen that there exists a critical point in the case of amplitude damping channel 
for initially unentangled state at which both firms have equal payoffs. The 
game transforms from leader advantage to the follower advantage game beyond 
this point, for highly decohered case in the presence of memory. For initially 
entangled state under the infiuence of amplitude damping channel we found two 
critical points. The game behaves as a follower advantage game within these two 
critical points. In the case of depolarizing channel the high correlation results 
in high payoffs. However, phase damping channel has no effect on the game 
dynamics. 
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2 Stackelberg duopoly game 



Stackelberg duopoly is a market game, which is rather different from the Cournot 
duopoly game. In Cournot duopoly game, two firms simultaneously provide a 
homogeneous product to the market and guess that what action the opponent 
will take. However Stackelberg duopoly is a dynamic model of duopoly game 
in which one firm, say firm A moves first and the other firm, say B, goes after. 
Before making its decision, firm B observes the move of firm A. This transforms 
the static nature of Cournot duopoly game to a dynamic one. Firm A is usually 
called the leader and firm B the follower, on this basis the game is also called 
the leader- follower model [30]. In classical Stackelberg duopoly it is assumed 
that firm B will respond optimally to the strategic decision of firm A. As firm A 
can precisely predict firm _B's strategic decision, firm A chooses its move in such 
a way that maximizes its own payoff. This informational asymmetry makes the 
Stackelberg duopoly as the first mover advantage game. 

A number of authors have proposed various quantization protocols for ob- 
serving the behavior of Stackelberg duopoly game in the quantum realm [9,29- 
32] . It has been shown that quantum entanglement affects payoff of the first 
mover and produces an equilibrium that corresponds to classical static form of 
the same game j^Sj. The effects of decoherence produced by various prototype 
quantum channels on quantum Stackelberg duopoly have been studied by Zhu 
et al. [20]. We study the effects of correlated noise on quantum Stackelberg 
duopoly, using amplitude damping, phase damping and depolarizing channels. 



3 Calculations 

In a quantum Stackelberg duopoly game, for each firm A and B the game space 
is a two dimensional complex Hilbert space of basis vectors |0) and that is, 
the game consists of two qubits, one for each firm. We consider that the initial 
state of the game is given by 

IV'i) = cos6i|GG) +sin6l|ll) (1) 

where is a measure of entanglement. The state is maximally entangled at 
= J. In the presence of noise the evolution of an arbitrary system can be 
described in terms of Kraus operators as dJ 

I 

where — \il>i){ipi\ is the initial density matrix and the Kraus operators Ei 
satisfy the following completeness relation 

Y^ElEi^l (3) 

I 
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The single qubit Kraus operators for uncorrelated quantum amplitude damping 
channel are given as [23l [H] 
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The Kraus operators for amplitude damping channel with correlated noise for 
a two qubit system are given as [24] 
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The action of such a channel on the the initial density matrix of the system is 
given by 

1 1 

P=(1-M) E El^P^E^n+^^Y.^uP^K^ (6) 
m,n— 

where the superscripts u and c represent the uncorrelated and correlated parts 
of the channel, respectively. The above relation means that with probability 
fj, the noise is correlated and with probability (1 — fi) it is uncorrelated. The 
Kraus operators for phase damping channel with uncorrelated noise for a system 
of two qubits are given as [23, 24 



= Veme„cr,„ (E) (Jn, m,n = 0, 3 

whereas the one with correlated noise are given as 



E?, 



eiCTi (K) (Ti, 



1 = 0, 3 



(7) 



(8) 



Similarly, the Kraus operators for depolarizing channel are described by equa- 
tions ([7]) and ([8]) with indices run from to 3. Where ctq is identity operator for 
a single qubit and cti, a2 and (73 are the Pauli spin operators. For phase damping 
channel, eo ~ (1 — pi) and 63 = pi and for depolarizing channel cq = (l—pi), 
ei = 62 = £3 = ^pi, where pi correspond to the decoherence parameters of the 
first and second use of the channel. The action of such a channel on the quan- 
tum system can be defined in a similar fashion as described earlier in equation 
®. 

In quantum Stackelberg duopoly game each firm has two possible strategies 
/, the identity operator and C, the inversion operator or Pauli's bit-fiip operator. 
Let X and 1 — x stand for the probabilities of / and C that firm A applies and 
y, 1 — y, are the probabilities that firm B applies, respectively. The final state 
after the action of the channel is given by |32j 



Pf 



= xylA®lBplA®ll+x{l--y)lA<?iCBplA 



'4 



+y{l-x)CA<ElB pC\®ll 

+ (1 - a;) (1 - y) Ca®Cb pC\® C 



(9) 
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where p (equation is the density matrix of the game after the channel action. 

Suppose that the player's moves in the quantum Stackelberg duopoly are 
given by probabilities lying in the range [0, 1]. In classical duopoly game the 
moves of firms A and B are given by quantities qi and 92, which have values 
in the range [0,oo). We assume that firms A and B agree on a function that 
uniquely defines a real positive number in the range (0, 1] for every quantity gi, 
q2 in [0, 00). Such a function is given by 1/(1 + qi), so that firms A and B find 
X and y, respectively, as 

^ y-T^ (10) 



1 + gi 1 + 92 

The payoffs of firms A and B are given by the following trace operations 

(qi, 92) = Tr[p^.P7 (91,92)], (gi, 92) =Tr[p^P°P (91,92)] (11) 

where P'^ , are payoff operators of the firms and are given by 

PA^{qi,q2) = — (fcPii - P22 - P33) 
912 

Pb^ (91, 92) = — (fcpii - P22 - P33) (12) 
912 

where are the diagonal elements of the final density matrix, fc is a constant 
as given in ref. [30] and 912 is given by 

^ (l + 9i)\l + 92) ^''^ 

The backward-induction outcome in the Stackelberg duopoly is found by 
first finding the reaction R2 (91) of firm B to an arbitrary quantity 91 chosen 
by firm A. It is found by differentiating firm _B's payoff with respect to 92, and 
maximizing the result for 91 and can be written as 

i?2 (91) = max Pb (91,92) (14) 

Once firm B chooses this quantity, firm A can compute its optimization problem 
by differentiating its own payoff with respect to 91 and then maximizing it to 
find the value 91 = 9i. Using the value of 9i in equation we can get the 
value of 92 . These quantities define the backward-induction outcome of quantum 
Stackelberg duopoly game and represent the subgame perfect Nash equilibrium 
point. The payoffs of the firms at this point can be found using equation (|lip . 



4 Results and discussion 

We suppose that before the firms measure their payoffs the game evolves twice 
through different quantum correlated channels. That is, the state, prior and 
after the firms apply their operators, is influenced by the correlated noisy chan- 
nels. 
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4.1 Correlated amplitude damping channel 

The subgame perfect Nash equihbrium point for the game under the influence 
of correlated quantum amplitude damping channel becomes 

* _ -fccOS^ 61 + ^1 {pi,p2,Hi,H2) 

-4 + A2 (pi,p2,Aii,Ai2) 

92 = -ifi , p / ^ (15) 

16 + B2[Pl,P2,Hl,H2) 

where the damping function Ai and Bi are given in appendix A. 

If we consider the influence of decoherence in the second evolution only 
{pi = 0), equation ([T5|) for unentangled initial state reduces to the following 
form 

1 

91 = 



-4^2(1-^^2) 

a* ^ l-2p2(l-M2) . . 

^2 2[2-p2{6-5A*2-P2(l-M2)(5-8Ai2)}] ^ ' 

Here we have taken k — 1. The firms' payoffs under this situation become 



- 16p2(l - A*2) 



" 8{2 -p2(6 - 5^2) +Pi(8A*2 + 5)(m2 - 1)} 

In the classical form of the duopoly the perfect game Nash equilibrium is a 
point, whereas in this case it is a function of both decoherence and memory 
parameters. It can be easily seen that the results of ref. ^D] are retrieved by 
setting ^2 = 0, and setting P2 = reproduces the results of classical game. The 
existence of Nash equilibrium requires that firms' moves (ql and (72) should have 
positive values. It can easily be checked using equation ()16p that in the absence 
of quantum memory, no Nash equilibrium exists for P2 > ^- The presence of 
quantum memory allows the existence of Nash equilibrium for the entire range of 
values of p2, when /ij > 0.85. To see the influence of decoherence and quantum 
memory on firms' payoffs at the subgame perfect Nash equilibrium, we plot 
the payoffs (equation (|17p ) in figure 1 as a function of decoherence parameter 
P2. In the figure the dotted and dashed-dotted lines represent firm A and firm 
B payoffs for unentangled initial state, respectively. The superscripts u and 
e of Pa [Pb) in the figure stand for unentangled and entangled initial states 
respectively. It can be seen from the figure that a critical point exists due to the 
presence of memory at which both firms are equally benefited. This situation 
has not been observed, in the absence of memory, for unentangled initial state of 
the game. That is, in the absence of quantum memory the game is always first 
mover advantage game. For highly decohered channel and < 1, a transition 
from first mover advantage into second mover advantage occurs in the game 
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Figure 1: The payoffs {PA{ql,q2) and PB{qi,q2)) are plotted at the subgame 
perfect Nash equihbrium point against the decoherence parameter p2 both for 
entangled and unentangled initial states under the action of amplitude damping 
channel, with values of other parameters k = 1 and pi ^ fii — 0. For maximally 
entangled initial state the memory parameter — 0.9, and for unentangled 
initial state 112 — 0.7 



behavior. It can also be seen from equation (fT7|) that for fully correlated and 
fully decohered channel, both firms are equally benefited and get a payoff equal 
to |. It can also be shown that for smaller values of decoherence, the game is 
always first mover advantage irrespective of the value of the quantum memory. 

For a maximally entangled initial state of the game, the subgame perfect 
Nash equilibrium point becomes 

l-3pi(-l + /i2)~P2(2-3/.t2) 

4-8p2(l-A*2) 

[[-1+P2{2 + 3p2(-1 + A*2)-3a*2}] 
x{l + 2p2(-l + A^2)}] 

[-7 + P2(28 - 26/i2) + 9pt(-l + M2)' + pi 
X (-22 + 46/X2 - 23^1) - 6pi(2 - 5^2 + 3^1)] 
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The firms' payoffs corresponding to these values of q* become 



Pb 



Pa 



{-1 + P2(2 - 3^2) + 3pi(-i + 

32(l + 2p2(-l + M2)) 

[-{1 + 2p2(-l + M2)}{-1 +P2(2 - 3/i2) 



(19) 



[8{-7 + p2(28 - 26/i2) + 9p|(-l + fi^ f 
+p2(-22 + 46/i2 - 23m1) - 6p3(2 - 5fi^ + 3/il)}] 



One can easily check that these results reduce to the results of refs. [501 133] by 
setting ^2 = and P2 — respectively. The payoffs of firms for the maximally 
entangled initial state are plotted as function of decoherence parameter p2 in 
figure 1 . The solid line represents firm A payoff and the dashed line represents 
firm B payoff for maximally entangled initial state. One can easily check that 
for a maximally entangled initial state, the presence of quantum memory makes 
the firms better off as compared to the uncorrelated case of the channel. Fur- 
thermore, quantum memory shifts the critical point to higher payoffs than the 
critical point for maximally entangled initial state when the channel is uncor- 
related. This can also be shown that for a given decoherence level, the firms 
become worse off and the game becomes a follower advantage as the channel 
becomes more correlated. The effect of entanglement in the initial state on the 
firms' payoffs is shown in figure 2. Here we have taken into account the two 
decohering correlated processes. It can be checked that for certain range of de- 
coherence parameters, in the absence of quantum memory, the moves of firms 
are negative. This means that no Nash equilibrium exists in the region of these 
values of decoherence parameters. The presence of quantum memory, however, 
resolves this problem. As can be seen from figure 2, the presence of quantum 
memory in the payoff functions results into two critical points corresponding to 
two different values of entanglement parameter. The game is a follower advan- 
tage in this range of entanglement angle. The leader firm is worse off in this 
range of values of 9 and a global minimum in payoffs occurs at 9 — ^. 

4.2 Correlated Depolarizing channel 

When the game evolves under the influence of correlated depolarizing channel, 
the moves of firms at subgame perfect Nash equilibrium point become 



where the damping parameters A[, A'2, B[ and B'2 are given in appendix A. For 
unentangled initial state, firms' payoffs at the subgame perfect Nash equilibrium 



Slfccos^ 6 + A[ (pi,P2,A^i,A^2) 
3244 + v4^(pi,P2,/^i,/^2) 

6561fccos^ 9 + B[ (pi,P2,A*i,A*2) 
52488 + B^(pi,P2, Ml, M2) 



(20) 
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Figure 2: The payoffs of the two firms are plotted at the subgame perfect 
Nash equiUbrium point against the entanglement angle 6 under the action of 
amplitude damping channel, with other parameters pi = P2 = 0.5, k = 1 and 
/"i = = 0-5. 

point become 

p ^ (3-2p2)'{l + 2p2(-l + M2)P 

^ 24{-3-6p2(-l + M2)+4pi(-l + M2)} 

p ^ (3 - 2p2)'{l + 2p2(-l + A^2)}'{-3 + 2p2(-3 + 2p2)(-l + /i^)} 
48[9 + p2(-3 + 2p2){10 + 2p2(-3 + 2p2)(-l + M2)' " 9^2}] 

(21) 

where we have considered decoherence only in the second evolution of the game 
and have set A; = 1. It can easily be checked that for uncorrelated channel, firms' 
moves for p2 > 0.5 become negative and hence no Nash equilibrium exits for 
P2 > 0.5. On the other hand, it can be easily shown that the Nash equilibrium 
exits for the entire range of values of decoherence parameter, when the channel 
is correlated. The effect of quantum memory on firms' payoffs for unentangled 
initial state under depolarizing channel is shown in figure 3. The dashed line 
represents firm A's payoff and dashcd-dotted line represnets firm B's payoff for 
unentangled initial state. The payoffs grow up with increasing value of memory 




Figure 3: The payoffs {PA{qi,q2) and PB{qi,q2)) are plotted at the subgame 
perfect Nash equihbrium point against the memory parameter both for en- 
tangled and unentangled initial states under the action of depolarizing channel, 
with the other parameters pi ^ fii ^ 0, k ^ 1 and p2 = 0.5 



parameter and the game remains the first mover advantage game. 

When the initial state is maximally entangled and only the second decohering 
process is taken into account, the moves and payoffs of the firms for fc = 1 

1 3 

2 + 4{-3 + 2p2(-3 + 2p2)(-l + M2)} 

{~3 -t- 2p2(-3 + 2p2)(-l + M2)}{-3 + 4p2(-3 + 2p2)(~l + A^a)} 
63 + 8p2(-3 + 2p2){-9 + 2p2(-3 + 2p2)(-l + M2)}(-1 + ("2) 

(22) 

{3-4p2(-3 + 2p2)(-l + /^2)}' 
96{-3 + 2p2(-3 -I- 2p2)(-l + M2)} 

{3 - 4p2(-3 + 2p2)(-l + M2)}'{-3 + 2p2(-3 + 2p2)(-l + M2)} 
24[63 + 8p2(-3 + 2p2){-9 + 2p2(-3 + 2p2)(-l + M2)}(-1 + M2)] 

(23) 

From equation ([22]) one can easily check that no Nash equilibrium exists for 
decoherence parameter p2 > 0.32 when the channel is uncorrelated. However, 



become 

I1 = 

ql = 

Pa = 

Pb = 
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the presence of quantum correlations ensure the presence of Nash equihbrium for 
the entire range of decoherence parameter. It can be checked that for a given 
value of memory parameter the payoffs decrease when plotted as function of 
decoherence parameter and the game is always first mover advantage game. The 
dependence of payoffs for the maximally entangled initial state on the memory 
parameter is shown in figure 3. The solid line in the figure represents firm ^'s 
payoff and the dotted line represents firm i?'s payoff for the maximally entangled 
initial state. It can be seen that the game is no-payoff game around = 0.25. 
For other values of memory parameter, the game is first mover advantage game 
and the firms are better off when the channel is fully correlated. 




Figure 4: The payoffs (P^((7j[', (72) and PB{qi,q2)) arc plotted at the subgame 
perfect Nash equilibrium point against the memory parameter both for en- 
tangled and unentangled initial states under the action of depolarizing channel, 
with the other parameters pi = = 0, = 1 and p2 = 0.5 

To analyze the effect of entanglement in the initial state of the game on the 

Nash equilibrium, we plot the firms' moves q^, (j| in figures 4 and 5, respectively. 
From figure 4, one can see that the move of the leader firm is positive for the 
whole range of entanglement angle. However, for memory parameters Hi = 
II2 < 0.5, the move (72 of the follower firm is negative for a particular range of 
values of entanglement parameter as can be seen from figure 5. The negative 
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Figure 5: The move of firm B at the subgame perfect Nash equihbrium point 
under the influence of depolarizing channel is plotted against the entanglement 
angle 6 with other parameters pi = p2 = 0.25, fc = 1. 



value of q2 shows that no Nash equilibrium of the game exists in this range of 
values of 0. On the other hand, the existence of Nash equilibrium is ensured for 
the whole range of values of the entanglement parameter for highly correlated 
channel. The payoffs of firms as function of entanglement parameter are plotted 
in figures 6 and 7. A local maximum exists at ^ = tt for smaller values of memory 
parameters which disappears for the values of memory parameters > 0.5. 

4.3 Correlated phase damping channel 

The diagonal elements of the final density matrix of the game when it evolves 
under the action of correlated phase damping channel are given by 

Pii = 312 (cos^ + qiq2 sin^ 6) 

P22 = 912 (92 cos^ 9 + qi sin^ 0) 

P33 = 912 {Qi cos^ 6 + q2 sin^ 6) 

P44 = 912(9192 cos^6» + sin^6») (24) 
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Figure 6: The payoff of firm A at the subgame perfect Nash equilibrium point 
under the influence of depolarizing channel is plotted against the entanglement 
angle 9 with other parameters pi — P2 = 0.25, k — 1. 



The payoffs of firms in quantum Stackelberg duopoly game depend only on the 
diagonal elements of the final density matrix as can be seen from equations (jlll 
and [T^ . It is clear from equation that the diagonal elements of the final 
density matrix are independent of the decoherence parameters as well as from 
the memory parameters. Therefore, the correlated phase damping channel does 
not effect the outcome of the game. 



5 Conclusions 

We study the influence of entanglement and correlated noise on the quantum 
Stackelberg duopoly game by considering the time correlated amplitude damp- 
ing, depolarizing and phase damping channels using Kraus operator formalism. 
We have shown that in different entangling regions the follower advantage can 
be enhanced or weakened due to the existence of the initial state entanglement 
influenced by different correlated noise channels. The problem of nonexistence 
of the subgame perfect Nash equilibrium in various regions due to the presence 
of decoherence is resolved by quantum memory. 
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Figure 7: The payoff of firm B at tlie subgame perfect Nash, equilibrium point 
under the influence of depolarizing channel is plotted against the entanglement 
angle 6 with the following values of the other parameters pi = p2 = 0.25, k = 1. 



It is shown that under the action of amplitude damping channel, for initially 
unentangled state, the presence of quantum memory results into a critical point 
at the subgame perfect Nash equilibri\.im. The firms are equally benefited at 
this point and the leader advantage vanishes. Beyond this critical point, the 
Nash equilibrium of the game gives higher payoff to the follower firm as a result 
of quantum memory, that is, the game converts from leader advantage to the 
follower advantage game. Quantum memory, in the case of amplitude damping 
channel, favors the follower flrm in a particular range of values of entanglement 
parameter at the subgame perfect Nash equilibrium. In this range of values of 
entanglement parameter, the leader firm get worse off and the follower firm is 
better off (see figure 2). It is also observed that quantum memory validates the 
existence of Nash equilibrium for the whole range of entanglement angle and 
decoherence parameter. It is also shown that quantum memory in the case of 
phase damping channel has no effect on the subgame perfect Nash equilibrium 
and thus does not change the outcome of the game. 

In the case of depolarizing channel quantum memory and entanglement in 
the initial state influences the flrms' payoffs at the subgame perfect Nash equi- 
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librium strongly in a way different from amplitude damping channel. It is seen 
that for p > 0.32, no Nash equilibrium exists in case of unentangled initial 
state. Whereas the presence of entanglement in the initial state extends the 
span of decoherence parameter p from — 0.32 to — 0.5 for the existence of 
subgame perfect Nash equilibrium in the absence of quantum memory. On the 
other hand, we have observed that in the presence of memory, the subgame 
perfect Nash equilibrium exists for the entire range of decoherence parameter in 
both of the situations (for entangled and unentangled initial states). Similarly, 
it has been shown that memory has a striking effect that there exists a Nash 
equilibrium of the game for the entire range of entanglement parameter as well. 
Whereas, in the absence of memory, the noisy environment limits the subgame 
perfect Nash equilibrium to exist in a particular range of entanglement angle. 
In addition, a local maximum in payoffs is observed for small values of quantum 
memory parameters /Xj^, fj,2. For highly correlated channel this local maximum 
disappears and the payoffs reduce to zero. Unlike amplitude damping channel, 
the correlated depolarizing channel does not give rise to a critical point at the 
subgame perfect Nash equilibrium and as a result, the game always remains a 
leader advantage game. 
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Appendix A 

The damping functions Ai in equation (jl5|) are given as 

Ai = i[-A:-4(-l+pi)(-l+p2)p2(-l + M2) 

-2kp2(p2 + -P2^J■2) - 2fcpi[A*i +P2{2 -_p2(-2 + 

X (-1 + M2) - 4m2 + Mi(-2 + 3;U2)}] 

-A:cos26l + 4(-l +pi)(-l +P2)P2(-1 + M2)cos26l 

+ 2kp2{p2 + M2 ^ J52M2) cos 20 + 2fcpi[yUi +P2{2 - J32(-2 + Ail) 

x(-l + ^2) - 4^2 + Ml (-2 + 3,^2)}] cos 261] - 2pi(-l + ^1) 

x[2+J3i[-2 + fc{-l+p2(2-3/i2)+J32(-l+M2)} 

-4p2(-l + ^^2) + 2P2(-1 + A*2)] + 4P2(-1 + M2) 

-2p2(_i + /,2)]sin'e (25) 
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A2 = -4(l + fc)p2(-l + Ai2)-4(-l+Pi)pi(-l + Mi) 
X [-2 + k{-l+p2{2 - 3^2) +pI{-1 + /i2)} 
-4p2(-l + ^^2) + 2p^(-l + M2)] sin^ e (26) 

The damping functions Bi in equation (jlSp are given as 

Bl - -2+ifc[l + 2p2(p2+/i2-P2M2) 

+2piK +P2{2 - 4^2 + Ail (-2 + 3^2)}] 

+ [1 - 2[P2(P2 + A*2 - ^2^*2) +Pi{Aii +P2(2 - 2^1 

-4^2 + 3AiiAi2)}]] cos 261] - \pi{~l + Mi)[fc(-4 + 3pi) 

x{-l +p2(2 - 3^*2) +P2(-1 + M2)} + 6(-l +pi){-l - 2p2 

X (-1 + ^l2)+vl{-l + ^^2)}] sin' e - \p2{-l + ^^2) 

X [4(1 + k) + [2(-l + P2) +pi{2 + p2(-2 + fc(-2 + ^1))}] sin2 9] 

(27) 

B2 = -16 + 2[2{-l-(l + fc)p2(-l + M2)}+P?(-l + Mi) 

X [(2 + fc)(-l + J52)' +P2{4 + 3fc - (2 + k)p2}^l2] 

+Pi{~l + /Xi)[-(2 + fc)(-l +P2)' +P2{-4 - ik 

+ (2 + k)p2}^i2] + {-l+pi)pi{-l + Mi)[-(2 + fc)(-l +P2)' 

+P2{-4 - 3fc + (2 + fc)p2}M2] cos 261] [{1 + (1 + k)p2 

X(-1 +^2)}C0S^6I+ [1 + (1 + /s)p2(-l +/i2) 

X [(2 + fc) (-1 + P2)' + P2{4 + 3fc - (2 + fc)p2}/i2] 
+pI{-1 + Mi)[-(2 + k){-l + p2f + P2{-4 - ik 
+ (2 + /c)p2}/i2]sin' 0] - [-2(-l+pi)[pi(-l + Mi) 

X{-1 + (-2 +P2)P2(-1 + M2)} + (-1 +P2)P2(-1 + M2)] 
-fc[l +P2(P2 +M2 -P2/i2) +Ml){-(-l +P2)' 

+ (-3 +P2)p2Ai2} +P2{2 - P2(-2 + /^i)(-l + ^^2) 

-A^i^ + /ii(-2 + 3^2)}]] + (-1 +Pi)[2{pi(-1 + Ml) 

X{-1 + (-2 +P2)j52(-1 + M2)} + (-1 +P2)j52(-1 + /i2)} 

+ fc[(-l +P2){-1 +P2(-1 +M2)} +/^l){-(-l +P2)' 

+ (-3 + P2)p2Ai2}]] COS 20][p2{-2 - (2 + fc)j52(-l + M2) 

+ (2 + /c)/i2} cos^ 61 + [2j9i[-l + +J3i(-1 + Ml) 

X{-1 + (-2 +P2)P2(-1 + M2)} - P2{1 + (-2 +P2)Mi}(-1 + 1^2)] 
+ fc[l +Pl[-2 + /ii +pi(-l + /ii){-(-l +P2)' + (-3 +P2)P2M2} 

+P2{2 - 2/^2 + /^i(-2 +P2 + 3m2 - P2A*2)}]]] sin" Q] (28) 
The damping functions A[ in equation ()20|1 arc given as 
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a; = 4(2 + fc)pi(-3 + 2pi)(-l + ^i){-9 + 8p2(-3 + 2p2) 
X (-1 + ^2)} - 36(2 + fc)p2(-3 + 2p2)(-l + M2) 

+ ^fc[9 + {9 - 24p2 + 8pi(-3 + 4p2)} cos 20] (29) 

A'2 = 8(2 + fc)[pi(-3 + 2pi)(-l + Aii){-9 + 8p2(-3 + 2p2) 

X (-1 + /i2)} - 9p2(-3 + 2p2)(-l + Ih)] (30) 

The damping funetions B[ in equation (j20p are given as 

- 8(2 + A:)2(3-2p2)'p2{9-8pi(-3 + 2pi)(-l + /ii)}2 
X [fc{-9 + 4pi(-3 + 2pi)(-l + Ml)} + 8pi(-3 + 2pi) 
X (-1 + Ml) - A:(-3 + 4pi)(-3 + 4^2) cos 26*] + 18(2 + k) 
xp2{-3 + 2p2){-9 + 8pi(-3 + 2pi)(-l + Aii)}(-1 + M2) 
X [(2 + k){9 - 8pi(-3 + 2pi)(-l + Ml)} + k{-3 + 4pi) 
X (-3 + 4^2) cos 20] (31) 

B!, = -72(2 + fc)p2(-3 + 2p2)[4{-9 + 2pi(-3 + 2pi)(-l+Mi)} 
+fc{9 + 4pi(-3 + 2pi)(-l + Mi)}]{-9 + 8pi 
X (-3 + 2pi)(-l + Mi)}(-1 + M2) + 16(2 + kf{3 - 2^2)' 
xpl{9 - 8pi(-3 + 2pi)(-l + Mi)}'(-1 + 
+81[fc2{-81 + 8pi(-3 + 2pi){9 + 2pi(-3 + 2pi)(-l + Mi)} 
x(-l + Ml)} + 64pi(-3 + 2pi){-9 +pi(-3 + 2pi) 
X (-1 + Mi)}(-1 + /ii) + 16fcpi(-3 + 2pi){-9 + 4pi 
X (-3 + 2pi)(-l + Mi)}(-1 + /^i) + fc'(3 - 4pi)2 
x(3-4392)^cos2 26i] (32) 
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